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Abstract 

We obtain a criterion for pulsating front speed-up by general periodic incompressible flows in 
two dimensions and in the presence of KPP nonlinearities. We achieve this by showing that the 
ratio of the minimal front speed and the effective diffusivity of the flow is bounded away from zero 
and infinity by constants independent of the flow. We also study speed-up of reaction-diffusion 
fronts by various examples of flows in two and three dimensions. 

1 Introduction 

We consider reaction-diffusion fronts propagating in a strong periodic incompressible flow on M. n : 

T t + Au-VT = AT + f(T). (1.1) 

Here T(t,x) G [0, 1] is the normalized temperature and the nonlinearity / is of the KPP type: f(s) 
is a Lipschitz function such that /(0) = /(l) = 0, while f(s) > and f(s) < f'(0)s for s G (0, 1). 
The Lipschitz flow u(x) is 1-periodic, incompressible, and has mean zero. That is, u(x + k) = u(x) 
when k G Z™, V • u = 0, and 

/ u(x)dx = 0. 

The parameter A G K is the amplitude of the advection, and we will mainly be interested in A S> 1. 

It has been proved in [3] that when u[x) is periodic, equation (jl.ip has pulsating front solutions 
of the form T(t,x) = U(x ■ e — ct,x), where c > is the propagation speed and e G M. n is the unit 
vector in the direction of propagation. The function U(s,x) is periodic in x G 1" and has uniform 
in x G T n limits as s —> ±oo: 

lim U(s,x) = 1, (1.2) 

s— 00 

lim U(s, x) = 0. 

Pulsating front solutions were shown in [3j to exist for all |e| = 1 and all c > c*(A). As in the 
one-dimensional case without advection [15], the minimal front speed c*(A) (we suppress the u and 
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/ dependence in our notation) determines the propagation speed of solutions of the Cauchy problem 
for (jl.ip with general compactly supported initial data, and is therefore of a special interest [H [18] . 

The presence of an incompressible flow in improves mixing due to diffusion and is thus 

expected to enhance the speed of reaction-diffusion fronts. This problem has been studied actively 
in the recent years, especially in the large A limit: it has been shown in [2J El Q2] that the pulsating 
front speed in the direction of a mean-zero shear flow behaves as c*(A) = 0{A) for large A, and in [T7] 
that c* e {A) = 0{A l / A ) for cellular flows in two dimensions. In both of these cases the minimal front 
speed scales as c* e {A) ~ \J D e (A) for A 3> 1 where D e (A) is the corresponding effective diffusivity of 
the flow Au in the direction e. In the present paper we use the method of [17] to show that this is 
indeed the case in general in two dimensions. 

Let us recall the definition of the effective diffusivity. Consider the advection-diffusion problem 

p t + Au ■ Vp = Ap (1.3) 

with u a periodic incompressible flow. The long-time behavior of the solutions of (jl.3p is governed 
by the effective diffusion equation 

i,3=l J 

The (x-independent) effective diffusivity matrix cr(A) is obtained as follows. For any e S R n , let 
Xe(x) be the periodic mean- zero solution of the cell problem 

- Axe + Au ■ Vx e = Au ■ e (1.5) 

on T n . Then the matrix a (A) is given by 

e-a(A)e'= [ (V X e + e) • (Vxe' + e')dx = e ■ e' + / V Xe • VXe'dx, (1.6) 

for any e, e' G M. n . The effective spreading in the direction e is then governed by the effective 
diffusivity 

D e (A) = e • a(A)e = 1+1 \Vxe\ 2 dx. 

When the nonlinearity in (jl.ip is weak and (jl.ip becomes 

T t + Au ■ VT = AT + ef(T) 
with e<1, one may consider the long time-large space scaling t — > t/e 2 , x — > x/e leading to 

T t + ju(^yVT = AT + f(T). 

The homogenized version of this equation is 

T t = V-(cr(A)VT)+f(T), 

with the corresponding homogenized minimal front speed v*(A) = 2W f'(0)D e (A). This approxi- 
mation holds only on certain time-space scales in the original variables (namely, t = 0(l/e 2 ) and 
x = 0{l/e)). However, it suggests a relation between the minimal front speed for (jl.ip and the 
effective diffusivity. The following result confirms this relation in two dimensions. 
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Theorem 1.1. There exists C > (independent of A,u, f,e) such that if u(x) is a 1-periodic 
incompressible Lipschitz flow onM 2 , f a KPP nonlinearity, A £ R ; and e£l 2 a unit vector, then 

< <(A) < C y/77W 1 + ). (i.7) 

C{l + ^fW)) ~ ^DjA) ~ VJKn VJy>> 



Moreover, there is fo > such that 



C* e (A) 



2y/f'(0)D e (A) 



<c/'(o) 



1/4 



whenever < /'(0) < /q. 



That is, the ratio c*^)/^ D e (A) is bounded away from zero and infinity by constants only 
dependent on /'(0), and becomes close to 1\J /'(0) when f'(0) is small. We note that the slightly 
weaker, and only upper bound c* e {A)/ ^ D e (A) < C eV / /'(0)(l + ■ s /f(fi))\\Au\\l for any e > has 
been obtained in [12j . 

As /'(0) — > +oo, the lower bound in (|1.7p stays bounded whereas the upper one grows linearly 
with /'(0). We show by looking at the example of shear flows (see Example 13. 2\i below that at least 
the lower bound cannot be improved. This conclusion can also be reached using the results of [12] 
for the shear flow. 

It follows from Theorem 1 1 . 1 1 that the minimal front speed c* e {A) has the same asymptotic behavior 
in the regime of large A as does ^D e (A). But for the latter quantity we have the following general 
sharp criterion which holds in any spatial dimension. 

Proposition 1.2. Let u(x) be a 1-periodic incompressible Lipschitz flow on M n and let e G M n be a 
unit vector. 

(i) If the equation 

has a solution 4> e G ff 1 (T n ) ; then 



u ■ V4> e = u-e (1.9) 



limsupL» e (yl) < +oo. (1.10) 

(ii) If hi. 9^ has no H 1 (T n ) -solutions, then 

lim D e (A) = +oo. (1.11) 

Remark. Note that it follows from part (i) that the set of all non-negative multiples of unit 
vectors e E W 1 for which (jl.lOj) holds is a subspace of M n . Indeed, the sum of two solutions of (|1.9|) 
for e and e' is a solution for e + e', and the negative of a solution for e is a solution for — e. 

The result of Proposition 11.21 is not new. It has already appeared in [llj, although it has been 
stated only in two dimensions, and the first claim has also appeared earlier in jJJJ. Much more precise 
asymptotic behavior of D e (A) is well understood for many specific examples of flows — see [16] for 
an extensive list of references. 

Putting Theorem 11.11 and Proposition 11.21 together, we have the following characterization of 
flows in two dimensions that speed up KPP fronts. 
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Corollary 1.3. Let u(x) be a 1-periodic incompressible Lipschitz flow on M?, let f be a KPP non- 
linearity, and e G R 2 a unit vector. 

(i) If CCH has a solution cf> e G H 1 ^ 2 ), then 

limsupc*(A) < +00. (1.12) 

(ii) If (|1.9p has no H 1 ^ 2 )- solutions, then 

lim c*(A) = +00. (1.13) 

Remark. In particular, the pulsating front speed for KPP nonlinearities may not diverge to +00 
along some sequence of amplitudes while staying bounded along another sequence. 

The paper is organized as follows. Section [2] contains the proofs of Theorem 11.11 and Propo- 
sition 11.21 The proof of Theorem 11.11 is based on the proof of the main result of [T7] . Section [3] 
contains the generalization of Corollary |1.3( ii) to higher dimensions and various examples. In par- 
ticular, we show there that the minimal front speed for a class of cellular three-dimensional flows 
satisfies (|1.13p . To the best of our knowledge, this is the first time that the front speed-up by a 
cellular three-dimensional flow has been established. 

Acknowledgment. This work has been supported by ASC Flash Center at the University of 
Chicago. LR was supported by NSF grant DMS-0604687 and AZ by NSF grant DMS-0632442. 

2 Diffusivity enhancement and front speed-up 

We first present the proof of Proposition 11.21 for the convenience of the reader. 
2.1 The proof of Proposition [1.21 

Let us assume that there exists a sequence A n — > +00 and a constant M > such that D e (A n ) < M 
for all n. It follows from (jl.6p that there exists a sequence of mean-zero functions Xn( x ) on T n which 
satisfy 

- Axn + A n u ■ Vxn = A n u ■ e (2.1) 

such that HVxnlll — M for all n. As the functions Xn are uniformly bounded in // 1 (T n ) there exists 
a subsequence Xn k which converges to a function x( x ) weakly in H l (Y n ) and strongly in L 2 (T n ) as 
k — > +00. We divide (12. ip by A nk and pass to the limit k — > +00 to obtain 

u-Vx = u-e (2.2) 

in the sense of distributions. Since x G -ff 1 (T n ), p.2p holds almost everywhere on T n . This proves 
Proposition ll.2f ii) . 

In order to prove Proposition II .2( i) let us assume that a mean-zero function (j) e G if 1 (T Tl ) satisfies 
(jl.9p and let Xe be the mean- zero solution of (jl.5p . Consider the function r\ = Xe — (f>e which satisfies 

-A(rj - 4> e ) + Au ■ V77 = 0. 
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Multiplying the last equation by r/ and integrating by parts, from incompressibility of u we obtain 

J \Vn\ 2 dx = - J V0 e • Vndx. 
It follows that ||V?7||2 < ||V0 e ||2, and so 

D e (A) = l+[ \Vxe\ 2 dx<l + 2 \Vr]\ 2 dx + 2 \V(f) e \ 2 dx < 1 + 4||V<£ e |||. 
Therefore, D e {A) is uniformly bounded in A and (i) follows. □ 
2.2 A variational principle for c*(A) 

The proof of Theorem 1 1 . 1 1 relies on a variational principle for the effective speed which we now recall. 
Details and proofs of Propositions 12.11 and 12.21 below can be found in [31 SI El H3 E] • Consider the 
eigenvalue problem on T n 

A<p — Au ■ Vip — 2Ae • Vp + XAu • eip = K e (X; A)tp, ip > 0. (2.3) 

It has a unique eigenvalue n e (\; A) that corresponds to a positive periodic eigenfunction (p e (x; A, ^4). 

Proposition 2.1. The minimal front speed is described by the variational principle 

w . H a±*!+!M. (2.4) 

A>o A 



It is convenient to rewrite the eigenvalue problem (J2T3J) in terms of the function 

^ e (x) = ip e (x)e Xx - e 

on M n , with ip e periodically extended from T n . This function is not periodic but rather belongs to 
the set 

E^ x = {ip(x) | ip(x)e~ Xx ' e is 1-periodic in x and ip > 0} . 
The corresponding eigenvalue problem for ip e { x ) is 

£V := - Au ■ Vip = fj, e (\; A)ip, ip € E+ x (2.5) 

with /j, e (X;A) = A 2 + K e (X;A) the unique eigenvalue of (|2.5p . The variational principle (|2.4p may 
now be restated as 

<(A) w/ilfM. (2.0) 
A>0 A 

Let us recall some basic properties of the function fX e (X; A). These can be found, for instance, in [31 
Proposition 5.7] (with our fj, e being their —h): 

Proposition 2.2. For each fixed A £ R we have ^ e (0;A) = 0, and the function jj, e {X;A) > is 
monotonically increasing and convex in A > 0. 
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Proposition 12.21 allows us to define the inverse function X e (fi;A) to (j, e (X;A) that is increasing 
and concave in \x > for a fixed A. The eigenvalue problem (|2.5p may be now re-formulated as 
follows: given \x > find A = A e (/i; A) so that the problem 

Atp - Au ■ = /x^, (2.7) 

has a solution ^ £ -^a - ^he variational principle ([2.60 for the minimal front speed now becomes 

c e A) = mf — — — . (2.8) 
m>o X e (p;A) 

2.3 The proof of Theorem [TTT1 

The proof of Theorem II. H is based on the ideas of [17] where they were used to obtain the asymptotics 
of the pulsating front speed for cellular flows in the direction (1,0). Here we extend them to general 
flows and directions e, and show that they yield the conclusion of the theorem. The main ingredient 
is the following Lemma. 

Lemma 2.3. There exists p > such that when e G [0, ^] and po{e) = pe 4 , then for any fx 6 
[0, //o(e)], any flow u as in Theorem any unit vector e £ M 2 , and any A £ R we have 

(1 - e)-£= < A e ( M ; A) < (1 + ^)^=. (2.9) 

We postpone the proof of this Lemma and prove Theorem 11.11 first. Fix any e > 0. Then the 
variational principle for c* e {A) implies that 

cJA) = mi -— — > mm<^ mf -— — , mf — — } (2.10) 

c>» A e (/i;4) lo<^<Me) X e (fi;A) ,j>^{e)X e {n;A)} 

and 

c* e (A)< inf {^^f . (2.11) 
In addition, as the function X e (fi,A) is concave and increasing in A and A e (0; A) = 0, we have 



At>Mo(e) A e (/i; A) A e (/i Q (e);A) l + £ 

It follows now from (|2TTUj) and (gSJ) that 



e^)>V^) min ( if /m±A^\ ( , 12) 

1 + e [0<At<W)(e) V/X J 

Using (|2.1ip and (|2.9[) we also arrive at 



C - M) <^M*L inf ffl±if. (2.13) 

l-£ 0< M <w(e) 
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If now f'(0) < |/x Q (i), we take e = ^ 1 (4/ / (0)) < | to obtain 

/'(0)+H 



inf _ 

0<M<Mo(e) 

and so we get from (|2.12p and (|2.13p that 



277^^5 < c:(j4) < 2^77(5) V^S. 



(2.14) 



Therefore, we have 



2 V //'(0)A ! (A) 



< 



1-e* 



for some C > 0. This proves the second claim in Theorem II .li 

On the other hand, if /'(0) > 3/^0(2)1 then we set e = 5 and obtain 



W/J l + l/2- eU - 1-1/2 
This and QHQ with s < \ for /'(0) < |/x (|) 

proves the first claim in Theorem I l.li □ 

2.4 The proof of Lemma 12.31 

The proof of Lemma 12.31 is similar to that of Theorem 3.1 in [T7j but our result is slightly sharper. 
First, we rewrite the eigenvalue problem (|2.7p in terms of the function (,{x\, x 2 ) = ln^(xi, 3:2): 



A(-Au- VC = /U- |VC| 2 , 

C(xi + 1, x 2 ) = C(x 1 ,x 2 ) + A e (/u; A)ei 

((xi,x 2 + 1) = ((xi,x 2 ) + A e (/u; yl)e 2 , 



(2.15) 



where e = (ei,e2). Without loss of generality we may assume that £(^1,^2) has mean zero on T 2 
(which is now viewed subset of M 2 ): 



Q{xi,x 2 )dx\dx 2 = 0. 



As VC is periodic both in x\ and x 2 , and u\{x\, x 2 )dx 2 = for each x\ and f u 2 (xi, x 2 )dxi = 
for each x 2 (because u is 1-periodic, mean- zero, and incompressible), by integrating (|2.15p over T 2 
we obtain 

|VC| 2 dxidz 2 - (2.16) 



The Poincare inequality then implies 



/ \C\ 2 dxxdx 2 < Cfi. 
Jx 2 



(2.17) 
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The function C( x ii x 2) may be decomposed as 
C(xi,x 2 ) = K{V,A) ix-e 



Xe(xi,X 2 ) + S(xi,X 2 ). 



1 



Here Xe(xi,x 2 ) is the mean-zero periodic solution of (|1.5p and S(xi,x 2 ) is a mean-zero periodic 
correction that we would like to show to be "small". Let us set 

1 



This function satisfies 



x 2 ) = A e (>; A) ( x • e - - - Xe(a?l, x 2 ) 



A$ - ylti • V$ = 

$(xi + 1, x 2 ) = x 2 ) + A e (/x; A)ei 
x 2 + 1) = *(a?i,x 2 ) + A e (^; A)e 2 . 



(2.18) 



The definitions of $ and D e (A), periodicity of Xe> and |e| = 1 imply that 



(2.19) 



\\V$h = \ e (fi;A)y/D e (A). 

Lemma 2.4. There exists a universal constant C > smc/i t/iai /or any flow u as in Theorem M.l 
any unit vector eGl 2 , and any A £ 1, we have 



1*1 



L°°(T 2 ) 



<CAe(M;^)-\/AP) 



(2.20) 



and 



C(^l,^2)<CVM for (xi,x 2 ) G T 2 . (2.21) 
We postpone the proof and first finish that of Lemma 12.31 It follows from (|2.20p and (|2,2ip that 

S(x 1 ,x 2 ) < c[x e (fi;A)y / D e (A)+^Jlj . (2.22) 

Note that this is only a bound from above, as is the one in (I2.2ip . On the other hand, as 

fi= / \X7(\ 2 dx 1 dx 2 = / \V(<$> + S)\ 2 dxxdx 2 , 
Jt 2 Jt 2 

the triangle inequality implies that 

IVM- ||V*|| 2 | < ||V5|| 2 . (2.23) 
It follows from (I2.15P and (I2.18P that the function S is a mean-zero periodic solution of 

AS - Au-VS = fi- |V($ + S)\ 2 . 
Multiplying both sides by S and integrating over T 2 we obtain using (|2.22l) . 

\VS\ 2 dx 1 dx 2 = J S|V($ + S)\ 2 d Xl dx 2 < C \x e (^A) y / D e (A) + ^Jl] [ |V(* + S)\ 2 d Xl dx 2 

= C 



X e ( y ^A) v / D e (A) + ^Jl 
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with C from Lemma 12.41 Using this inequality, (|2.23j) , (|2.19p , and \/ab + c < 5a + 5 1 b + ^/c with 
5 = e/2C, we obtain for any e > 

y/Ji < X e (fi; A)y/DjA) + Cy/\ e (n; A)^/DjA)fi + < (l + ^j A e (/z; Aj^+^p + C^ 4 . 
We now let /x (e) = (e/4C) 4 for e G [0, ±] so that 

^ + cy/ 4 < ^ 

fc o 

for all fi G [0, /J-o(e)]. Then for /i G [0, /Uo(e)] we have 



yfil<(l+e)\e(KA)y/D e (A). 

In a similar manner, from 



V£ > A e (/i; A) VAP) - Cy A e (/i; Ajv^pj/i + /i 3 / 2 
we obtain for G [0, ^o(e)], 

This finishes the proof of Lemma 12.31 □ 
2.5 The proof of Lemma 12.41 

Both statements of this lemma are an immediate consequence of the following two propositions. 

Proposition 2.5. There exists a universal constant C > such that if a G M 2 and the function q 
on M? satisfies 

Aq-Au- Vq+\Vq\ 2 > 0, 
qfa + 1, x 2 ) = q(xi,x 2 ) + a%, 
q{x x ,x 2 + 1) = q(xi,x 2 ) + a 2 , 

/ g(xi,x 2 )dxidx 2 = 0, 

then 

q(xi,x 2 ) < C(||Vg|| L 2 (T2) + \ ai \ + \a 2 \) (2.24) 

for all (xi,x 2 ) G T 2 . 

Remark. When a 2 = then this is Proposition 4 in [1TJ . 
Proof. The Poincare inequality says that 

IMIzW < (§-l)l|Vg|| L2(T2) (2.25) 
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for some C > 2. We will show that (|2.24p holds with this C. 
Let 

M : = max q(x 1 , x 2 ) = q(x ) 

(xi,x 2 )GT 2 

for some xo G T 2 and consider the case 

M > 2||g|| i 2( T 2) + |«i| + |cn 2 | 

(otherwise we are done by (|2.25p ). 
Assume first that 

for any x\ G [0, 1] there exists s(x±) G [0, 1] such that q{x\, s(x\)) > M — \a%\ — \a 2 \ 
Observe that if we define the set A C [0, 1] as 

A= {x 1 £ [0, 1] : 3r(xi) G [0, 1] such that \q(xt,r(xi))\ < 2||g|| i 2 (T 2)} 

then 



\A\ > \. 



(2.26) 

(2.27) 
(2.28) 



Now, (p^6]) - ([238l) imply that 



f 


dq 


2 

dX2 > 


[ -p-dx 2 


h 


dx 2 




.JO OX 2 



for all x\ G A. We then obtain 
|V 'q\ 2 dx\dx2 > 





dq 


L 


dx 2 



> (M -\ai\- \a 2 \ - 2||g|| L 2( T 2)) 



1 2 

dx\dx2 > - \M - \a\\ - \a 2 \ - 2||g|| i 2( T 2)) . 

It follows from the above and (|2.25p that 

max q{xi,x 2 ) = M < 2 (|| Vg|| L 2 (x2) + ||g|| L 2 (T2) + |ai| + \a 2 \) < C (||Vg-|| L 2 (T 2) + \a x \ + \a 2 \) ■ 

(xi,X2)GT J 

If ([236]) does not hold, but 

for any X2 G [0, 1] there exists s(x2) G [0, 1] such that q(s(x2), X2) > M — \a\\ — \a2\1 (2.29) 

then an identical argument applies. Assume therefore that neither (|2.26p nor (|2.29p hold. That 
means that there are £i>£2 G [0, 1] such that for all ^ G [0, 1], 

<M-\ ai \ - \a 2 \. 

If now Sj, j = 1, 2, 3, 4, are the four unit squares with sides parallel to the axes, and having (£1, £2) 
as a vertex, then for each j we must have 

max q(x\, x 2 ) < M — \a\\ — |a 2 | + |ai| + |«2| = M 

(xi,x 2 )G95j 

by the assumptions on q. But one of these squares contains xo with g(xo) = M, which is a con- 
tradiction because q satisfies the maximum principle. Hence (|2.26p or (|2.29p hold and we are done. 

□ 
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Proposition 2.6. For all A we have 

Xe(jji;A)<y/ii. (2.30) 



Proof. This is an immediate consequence of (12,15p . (12,16p . and |e| = 1: 

M = / iCrJ 2 + \( X2 \ 2 dxidx 2 > ( / Ca;i^l^2 ) +( / Cx 2 dxidx 2 ) = A e (/j; A) 2 . □ 



Proof of Lemma 12.41 This holds with C being (1 + \/2)-times the constant from Proposition! 
The upper bound on &(xi,x 2 ) follows directly from (|2.18p . (|2.19p . and D e (A) > 1. The lower bound 
on <3? follows after applying Proposition 12.51 to — $. The upper bound (I2.2ip for C,{ x i-, x 2) follows 
from Propositions 12.51 and 12.61 and (|2.16p . □ 

3 An extension to higher dimensions and examples 
3.1 The main result in higher dimensions 

As mentioned before, Proposition 11.21 holds in any dimension. Moreover, the following partial analog 
of Corollary 11.31 holds. 

Theorem 3.1. Let u(x) be a 1-periodic incompressible Lipschitz flow on R n , let f be a KPP non- 
linearity, and e 6 W 1 a unit vector. If (II. 9p has no H 1 (T n ) -solutions, then 

lim c*JAu,f)= lim D e (Au, f) = +oo. (3.1) 

Proof. Assume that there is a sequence A k — * +oo such that 

M := sup c* e {A k ) < +oo. (3.2) 

n 

Let <p k be a solution of (|2.3p on T n with A = A k and A = X k where X k is such that 

<{At) - m±*±^iiM ( 3. 3) 

Afc 

with K e defined in (|2.3p . Moreover, we choose (p k so that u) k = ln<£>& is mean-zero on T n . Denote 
K k = K e (X k ;A k ) so that 

Aw fc + \Vuj k \ 2 - A k u ■ Vuj k - 2X k e ■ Vu k + A k \ k u ■ e = n k . (3.4) 

Integrate this over T n and use that u is incompressible and mean-zero to obtain 

/ \Vu k \ 2 dx>0. 



K k 



This, (|3.2p . and (|3.3p mean that the A& are bounded away from and oo (namely, |2Afc — M\ < 
yj M 2 — 4/'(0) ). So after passing to a subsequence we can assume X k — ► Ao £ (0, oo). This and (|3.2p 
give 

sup ||Vo; fc ||^2(- T n) =supK fc < +00. 

fe K ' k 



11 



By (|2.25p . the oju are bounded in H (T n ), and so converge to some ujq G H 1 (T n ) strongly in L 2 (T n ) 
and weakly in i? 1 (T n ). Then also Au>k — > A^o in the sense of distributions and |Vwfc| 2 are bounded 
in L 1 (T n ). Divide (|3.4p by j4j, and pass to the limit so that in the sense of distributions, 

— u ■ Vu>o + Xou ■ e = 0. 

Since ujq G i7 1 (T n ), we now see that (f) e = Uq/Xq & i? 1 (T n ) solves (|1.9p . yielding a contradiction. □ 



3.2 Examples 

We conclude with several examples of flows in two and three dimensions which speed up reaction- 
diffusion fronts. We will consider the case e = e\ = (1,0) (resp. e = e± = (1,0,0) in three 
dimensions) so that (|1.9p becomes 

u • V(p = ui, (3-5) 
but our analysis easily extends to other directions of front propagation. 

Example 3.2. Shear flows. These are flows of the form u(xx,X2) = (0(^2), 0) with mean-zero a. 
In this case (12.30 with e = e\ becomes 



Aip - (Aa(x 2 ) + 2X)ip xl + XAa(x 2 )(f = k(X; A)ip, ip > 0. 
The unique solution to this equation is of the form ip(x±,x 2 ) = ip(x 2 ) and satisfies 

¥x 2 x2 + XAa(x 2 )p = k(X; A)(p, ip > 0. (3.6) 
Obviously then k(X; A) = k(XA; 1) and so we have 



lim <a _ lim inf m±4^i) _ Um inf no _ inf a±*i). 

A^oo A A^ooA>0 XA A^ooA>0 A A>0 A 

Notice that this shows that cl(A)/A is decreasing to a positive limit. This has been established in [2] 
and an alternative (variational) characterization of the limit has been provided in [15]. Multiplying 
(|3.6p by ip and integrating over T 2 one obtains k(X; 1) < A||a+||oo with a+ the positive part of a. 
Since the operator on the RHS of (|3.6p is self-adjoint, it is easy to see that in fact lim^oo k(X; 1)/A = 
|| ck_|_ || oo . From this it follows that for each /'(0), 

r <( A ) ^ II II V v c e( A ) 

hm — — - < a+ 00 = hm hm 

A^oo A /'(0)-+oo A-^oo A 

On the other hand, (11.51) becomes 



-Xx 2 x2 = Aa(x 2 ) 

with x{ x li x 2) = x{ x 2), and it follows that D e (A) = l + ^4 2 7 2 with 7 = ||V X2 (— A X2 ) _1 a||2. We thus 
have 

Hm < IKIIoc 



^/D e {A) 7 

for any /'(0) and a ^ 0. This shows that the lower bound in (jl.7p is optimal up to a constant for 
large /'(0). Note that this example extends without change to any dimension. 
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If a solution 4> G H l (Y 2 ) of (|3.5p exists, then the function ^(^1,^2) = 0({^i},^2) — ^1 belongs 
to i? 1 (2T x T, 2T) (here 2T = [0,2] with and 2 identified, and {x} is the fractional part of x\ in 
three dimensions ip G if 1 (2T x T 2 ,2T)). Moreover, ip satisfies 

u ■ Vf/> = and ^(^l + L ^2) = 4>{x\, £2) — 1- (3-7) 
In the following examples it will be a bit easier to work with t/j than with </>. 

Example 3.3. Percolating flows. Let u be such that for some X2 the streamline S of the flow 
(i.e., solution of the ODE X' = u(X)) starting at (0,2:2) reaches (1,^2) in finite time. This means 
that 5 is a periodic curve passing through (n, X2) for each nEZ, Such u is called a percolating flow 
(in the horizontal direction). Note that any non-zero shear flow is percolating. Since u 7^ on S, it 
is easy to show that if z/> from (|3.7p exists, it must be continuous and constant on S (see \19\ Lemma 
5.2]). This, however, contradicts the second condition in (|3.Tj) . Hence such a ip does not exist and 
Theorem 11.11 gives 

lim c* e (A) = +00 (3.8) 

for e = e\. This has been proved under some additional "non-degeneracy" conditions on u in [9]. 

Example 3.4. Cellular flows. These are flows with a periodic array of cells, each streamline of 
the flow being either a closed curve contained inside a cell or a part of the boundary of a cell. A 
typical example is the flow 

u(x 1 ,x 2 ) =V ± H(xi,x 2 ) = (-H X2 (x 1 ,x 2 ), H Xl (x 1 ,x 2 )) (3.9) 

with H(x\,X2) = sin 2irx\ sin 2irx2 the stream function of u. The streamlines of this flow are depicted 
in Figure [37TI 




Figure 3.1: A 2D cellular flow. 



13 



Constancy and continuity of ip from (|3.7p on each non-trivial streamline implies that ip has to 
be constant on the boundary of each cell (namely, it equals there the limit of the values on the 
streamlines approaching the boundary). But then ip can belong to H l only if it is constant on the 
whole "skeleton" of separatrices separating the cells, which again contradicts the second condition 
in (|3,7|) . Hence (|3.8[) holds for e = e\. We note that front speed-up by (certain generic) cellular 
flows has first been proved in [14J, with precise asymptotics established in [17] . 

Example 3.5. Checkerboard flows. Consider the cellular flows from the previous example with 
the flow vanishing in every other cell, thus forming the checkerboard pattern depicted in Figure [3721 




Figure 3.2: A checkerboard cellular flow. 

The above u is not Lipschitz in this case, so let us remedy this problem by taking, for instance, 
H(xi,X2) = (sin 2irxi sin 2irx2) a with a > 2 in the cells where u does not vanish. Now u vanishes on 
the boundaries of all cells, but if ip from (|3.T|) exists, the requirement ip £ H 1 still ensures that ip is 
constant on the boundaries of those cells in which the flow does not vanish (and it is continuous at 
these boundaries from inside of these cells). The values of these constants for two cells that touch 
by a corner must be the same. This follows from the fact that if an H l function has values a and 
b on two curves connecting a point x £ M 2 to a circle {y : \y — x\ = e} for some e > (and it is 
continuous at these curves, at least from one side), then a = b (see, e.g. |19[ Lemma 5.2] for this 
simple argument). We again have contradiction with the second condition in (j3.7l) . and (13. 8p for 
e = ei follows. 

Effective diffusivity enhancement for the flow in Figure 13.21 has been proved in [TT]. Notice, 
however, that our argument above can handle more general checkerboard-type flows. Consider, for 
instance, the case of such flows but with the contact between cells in which the flow does not vanish 
formed only by two touching cusps rather than two right angles. In this case the angle of contact 
between the cells can be equal to ir but we still have traveling front speed enhancement by the flow 
in the sense of (13. 8j) . 
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Example 3.6. Flows with gaps. Consider the cellular flow from Example 13.41 but with a vertical 
"gap" of width 5 > (in which u = 0) inserted in place of each vertical line {n} xl,n£Z, such as 
shown in Figure 13.31 




Figure 3.3: A cellular flow with gaps. 



This can be achieved by letting the stream function be, for instance, 

J (sin y^x\ sin27TX2) Q x\ mod 1 € [0, 1 — 5) 
jo xi mod l£ [1-5,1) 



H(x 1 ,x 2 ) 



with a > 2. More generally, we can assume that u has a periodic array of vertical channels in which 
the flow only moves "along" each channel. Namely, we assume that there is a C 1 map 7 : [0, 1] X R — » 
[0, l]xK which is one-to-one and satisfies 7(0, 6+1) = 7(0, b) + (0, 1) and u{^{a, b)) = c(a, 6)^7(0, b) 
for all (a, b) G [0, 1] x R and some c(a, b) £ R. This means that 7({a} x R) is a "vertical" streamline 
of u for each a (or a union of streamlines if u vanishes somewhere on this curve) . Let us also assume 
that the curve 7({0} x R) lies to the "left" of 7({1} x R). We now let u : [0, 1] [0, 1] be a C 1 
function with u>(0) = 1 and w(l) = 0. We define ^(7(0,6)) = uj{a) and set ip equal 1 and on the 
"left" and "right" components of [0, 1] xR \ 7([0, 1] x R), respectively. We have that ip is constant 
on each 7({a} x R) and so u • Vif) = 0. Since the periodicity of 7 ensures ip(xi,0) = ?/>(zi, 1), we 
can extend ip to an // 1 (2T x T, 2T)-function by letting ip{x\ + 1, x 2 ) = ip(xi,X2) — 1. Hence (j3.7|) is 
satisfied and (13.81) fails for e = e\. 



Note that the flows in the last example are percolating in the vertical direction if u does not 
vanish on each of the curves j({a} x R). Hence by the remark after Proposition II .21 and Example l3.3t 
speed-up of fronts in the sense of (|3.8p happens precisely when e ^ e\. On the other hand, (|3,8p is 
valid for any e in the case of the flows from Examples 13.41 and 13.51 

Our final example deals with speed-up of fronts by three-dimensional cellular flows. To the best 
of our knowledge this is the first example of this kind. 
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Example 3.7. 3D cellular flows. We consider here flows that have a cellular structure and are 
truly three-dimensional, with all three components of the velocity depending on all three coordinates. 
Such incompressible flows have been constructed in [7J. They have the form 

u(x 1 ,x 2 ,x 3 ) = (<5> xl (x l ,X2)W'{x 3 ),<S> X2 (x l ,X2)W'{x 3 ),k<S>(x l ,X2)W(x 3 )) (3.10) 

with A<I> = — k&. We will concentrate here on the simplest example of a flow with cubic cells. This 
flow is given by 

$(xi,X2) = cosxi COSX2, W(x 3 ) = sin x 3 (3-H) 

and k = 2, with W possibly any other 27r-periodic function vanishing at 0. 

The cube C = [0,2vr] 3 = (2vrT) 3 is a cell of periodicity for the flow (IXTTl) . Each of the eight 
dyadic sub-cubes of C of side-length tt is invariant under the flow and the flow in each of them has 
the same streamlines but moves along them in opposite directions in two neighboring sub-cubes. 
For the moment we restrict our attention to only one of them, [0, 7r] 3 . The streamlines of the flow in 
the part of this sub-cube given by x\ + x 2 > tt are depicted in Figure 13.41 (and they are symmetric 
across the plane x\ + x 2 = tt). 




Figure 3.4: A 3D cellular flow. 

Since 

ui(xi,x 2 ,X3) _ & Xl (xi,x 2 ) _ sinxicosx 2 
u 2 (xi,x 2 ,x 3 ) ® X2 (xi,x 2 ) cosxi sinx 2 

is independent of x 3 , the projection onto the (x±, X2)-plane of any streamline of the flow stays on a 
curve satisfying X' = V$(X). Figure [331 shows the phase portrait of this ODE in [0, tt] 2 . 

With the exception of the diagonals D = {(x±, x 2 ) : \xi — n/2\ = \x 2 — tt/2\}, each of the pictured 
curves connects a minimum and a maximum of $ (we will call these ^-curves). Any streamline of 
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Figure 3.5: A 2D projection of the 3D cellular flow. 

u in [0,7r] 3 , with the exception of those lying in (D x [0, 7r]) U ([0,7r] 2 x {0,7r}) (these form a set 
of measure zero and will be disregarded from now on), is a closed orbit, whose projection onto the 
(x\, x?)— plane is a portion of a $-curve contained in [0, tt] 2 \ D (in particular, it does not contain 
the corners of [0, 7r] 2 ). Some of these are stationary, namely those at (xi, x%, x%) with X3 = tt/2 and 
either x\ = ir/2 or X2 = tv/2. 

For each [x\,X2) G [0,vr] 2 \ D, let J(x\,X2) = (xi,ir — X2) if \x\ — vr/2| > \x2 — n/2\ and 
J(x\,X2) = (tt — x\,X2) if \x\ — tt/2\ < \x2 — tt/2|, so that (x\,X2) and J(x\,X2) lie on the same 
3>-curve. Extend J to the whole cube [0,7r] 3 by letting J(x\, X2, X3) = ( J(x\, X2), x-$) whenever 
\x\ — 7r/2| 7^ \x2 — 7r/2|. Notice that J is defined for almost all points in the cube (and we will from 
now on disregard the rest) and J 2 = Id. Each of the streamlines of u is symmetric either across the 
plane x\ = ir/2 or across X2 = vr/2, and hence the points x and J(x) lie on the same streamline of u. 

Let us now assume ip G i? 1 (2C,4vrT) with 2C = 4vrT x (2vrT) 2 satisfies 

u ■ = and ^(^1 + 2tt, X2, X3) = ijj{x\, X2, X3) — 2ir. (3.12) 

(which is the analogue of (|3.7p ) and assume without loss of generality that ip is real. Extend J to 
2C by J{x\ +ptt, X2 + qir, X3 + rir) = J(x\,X2,x^) + ir(p, q, r) for any p G {0, 1, 2, 3} and q, r G {0, 1} 
so that one still has that x and J(x) lie on the same streamline of u. The condition u ■ Vip = 
implies that tp is constant on almost all streamlines of u. In particular, for almost all x we have 
ip(x) = ip(J(x)). At the same time we know that the restriction of ip to almost every plane X3 = C is 
an H 1 (AttT x 2-zrT) function. This means that for almost every plane x% = C G (0, n) the restriction 
of ip to it (which we call VO belongs to i7 1 (47rT x 2itT) and satisfies ^(x) = ijj(J(x)) for almost all 
x = (xi,x 2 ) G 4vrT x 2vrT. 

Next we choose tq such that if = B((j<ir, 0), ro), then J B |V^(x)| 2 <ix < 7r/32 for jf = 0,1,2. 
It is then easy to show that for each j = 0,1,2 the set Rj of all r G (0, ro) such that 

- TT 

ess sup ip(jiT + r cos 9, r sin 0) — ess inf ip(j7r + r cos 0, r sin 9) < — 
e 2 

satisfies \Rj\ > 3ro/4. This is because 

— > f \Vii\ 2 dx > f [ 
32 Jb 3 J(0,r )\Rj Jo 



Ichp 
r 89 



d9dr > 



(0.n)\R, 



1 

2irr 



f-2-K 


dtp 


'0 


d9 



de 



> 



vr|(0,r )\^ 
8r 
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Let Ro = P\j =0 Rj so that |i?o| > r o/4, and let r G Ro. Then the values ip(x) for almost all points 
x on the circle Cj(r) = {(jvr + r cos 0, r sin 9) \ 9 G [0, 2-7r)} lie within an interval < vr/2. Since 

■0(r cos r sin #) = ^(7r — r cos r sin 0) and ^(7r + r cos r sin 9) = vp(2ir — r cos 9, r sin 9) for almost 
all (r, 9) £ -Ro x (— vr/4, vr/4) (because t/>(x) = '(/'(J(x))), we have that for almost all r G i?o> the values 
%j){x) for almost all x G U|=o^'j( r ) ne w ithin the interval Jo( r ) = U^=o^?'( r ) w ith |-^o( r )l ^ 37r/2. 
But this contradicts the second condition in (|3. 12|) because C2(r) = C${r) + (2tt, 0). Therefore there 
is no tp G H 1 {2C, 4ttT) which satisfies (|3.12p . and hence (|3.8p holds for e = e\. 

We note that our analysis can also be performed on other 3D flows from [7], for example, on the 
flow given by 

$(xi,X2) = 2cosa/3xi cos xi + cos ^2x2 — 7^ , W(xs) = sinx3, 

whose cells form a hexagonal 3D honeycomb lattice. Using the fact that $(xi,X2) = — <3?(i?(xi, X2)), 
with R the reflection accross any of the lines X2 = kx\ — ^, k = ±\/3, 0, one can show as above 
that the streamlines of the flow are symmetric accross the planes given by these three equations, 
and again conclude (|3,8p . 
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